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We study transitions between T states with the emission of charged pions using 604.6 fb" 1 
of data collected with the Belle detector at the KEKB asymmetric-energy e + e~ collider. The 
measured product branching fraction is B(T(iS) — > T(1S)tv + tv~) x B(T(1S) — ► (i + fi~) = 
(2.11 ± 0.30(stat.) ± 0.14(sys.)) x 10~ 6 and the partial decay width is r(T(4S) -> T(lS> + = 
(1.75 ± 0.25(stat.) ± 0.24(sys.)) keV. 

PACS numbers: 13.25.Gv, 14.60.Ef, 14.40.Aq 



The bottomonium state T(4S) has a mass above the 
threshold for BB pair production and decays mainly into 
5-meson pairs (B(?(4S) -► BB) > 96% Q). Recently, 
the decay modes T(45) — ► Y(mS)inr with m = 1,2 
as well as T(45) — > rjY(lS) j2] have been also ob- 
served. These decays as well as the anomalously large 
width of the T(55) -> T(mS)mr (m = 1,2,3) tran- 
sitions discovered by Belle Q give additional informa- 
tion about various QCD models that are used to de- 
scribe hadronic transitions of heavy quarkonia Q . Pre- 
liminary evidence for the decay T(45) — > T(lS')7r + 7r _ 
was presented by the Belle Collaboration in Ref. 
The BABAR Collaboration reported measurements of 
the T(45) transition to the T(15) or T(25) with the 
emission of a 7r + 7r~ pair 0, @. The first Belle mea- 
surements of T(4S) — > Y(lS)ir + ir~ were published in 
Ref. Q- The product branching fraction B(T(4S) — > 
T(15)7r+7r-) x B(T(1S) -» m + aO from Belle differs by 
2.4 standard deviations from the correponding value from 
BABAR In this paper we present a new study of the 
decay mode T(45) — > T(1S)tt + tt~ from the Belle ex- 
periment using a larger data sample and relaxed signal 
selection criteria. 

We use 604.6 ftr 1 of data collected on the T(4S*) 
resonance with the Belle detector Q at the KEKB 
asymmetric-energy e + e~ collider Q. We study Y(45) — ► 
Y(lS)ir + ir~ decays with a subsequent T(15) — > /i + £t~ 
transition. Charged particles arc reconstructed and iden- 
tified in the Belle detector, which consists of a silicon 
vertex detector (SVD), central drift chamber (CDC), 
aerogel threshold Cherenkov counters (ACC), time-of- 
fiight (TOF) scintillation counters, an electromagnetic 
calorimeter (ECL), and a A'i-muon detector (KLM). 

Charged tracks must originate from within a region 



of radius 1 cm and axial length ±5 cm centered in the 
e + e~ interaction point and not be associated with a well- 
reconstructed Kg meson, A baryon, or converted pho- 
ton; each charged track should have a momentum trans- 
verse to the beam axis (pr) of greater than 0.1 GeV/c. 
Charged particles are assigned a likelihood Li [13] (i = 
fj,, 7T, K) based on the matching of hits in the KLM to 
the track extrapolated from the CDC, and identified as 
muons if the likelihood ratio P^ = L^j^C^ + C n + Ck) 
exceeds 0.8, corresponding to a muon detection effi- 
ciency of approximately 91.5% over the polar angle range 
20° < 9 < 155° and the momentum range 0.7 GeV/c < 
p < 3.0 GeV/c in the laboratory frame. Electron iden- 
tification uses a similar likelihood ratio P e [ll| based on 
CDC, ACC, and ECL information. Charged particles 
that are not identified as muons and have a likelihood 
ratio P e <0.1 are treated as pions. Calorimeter clus- 
ters not associated with reconstructed charged tracks and 
with energies greater than 50 MeV are classified as pho- 
ton candidates. 



Candidates for T(45) 



the subsequent leptonic decay T(15) — > 



T(lS')7r + 7r _ decays with 
are se- 
lected from the standard hadronic-event sample for the 
first 492fb _1 data set (sample I), while an additional r- 
enriched sample is also used for the remaining 113 fb -1 
data set (sample II). The relevant selection criteria for 
the standard hadronic-event sample are the following: 
three or more charged tracks; a visible energy i? v j s of 
at least 0.2-y/s, where -y/i is the center-of-mass (cm.) en- 
ergy; a calorimeter energy deposit in the range 0.1 < 
E auxa /\fs < 0.8; and a maximum of 0.5^ for the mag- 
nitude of P z , the sum of the z components of the mo- 
menta of each charged track and neutral cluster, where 
the z axis is defined to be the direction opposite to the 
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positron beam. The variables E v [ s , E suul , and P z are 
evaluated in the cm. system. 

Unfortunately, the hadronic event sample criterion 
Esnm/y/s > 0.1 rejects a considerable fraction of our 
signal events. However, some of these lost events are re- 
covered in the T-enriched sample because it has a loose 
constraint on the -E sum variable (E sum < lOGeV). Only 
the newest r-enriched data are used because the ear- 
lier subset included a requirement on the sum of the 
magnitudes of the charged track momenta in the cm. 
frame (below lOGeV/c) that rejected most of our signal. 
Other selection criteria for the T-enriched sample, while 
not critical for this study, are enumerated below. The 
number of charged tracks in an event should be greater 
than one and less than nine with zero net charge. The 
maximum px among the tracks is required to be greater 
than 0.5 GeV/c. Beam-related background is rejected 
by requiring that the position of the reconstructed event 
vertex be less than 1 cm from the interaction point (IP) 
in the transverse direction and less than 3 cm from the IP 
along the beam direction. To suppress background from 
Bhabha and /i + (i~ events, the maximum opening angle 
between charge tracks is required to be less than 175° in 
the CM frame. 

To select T(4S) — > Y(IS)it + it~ decays, an event is 
required to contain exactly four charged tracks with 
a pair having an invariant mass M^ above 

9.0GeV/c 2 and a 7r + 7r~ pair whose opening angle V7r 
in the laboratory frame satisfies cos 6^ < 0.95. The lat- 
ter criterion suppresses the radiative return process [12 1 



T(15)7 as well as e H 



\i + \i 



7 wherein the 



photon converts to an e + e~ pair that is misidentified as 
a pion pair. Poorly reconstructed events are discarded 
by requiring a visible energy in the laboratory frame of 
10.5 GcV < < 12.5 GeV. 

To identify parent resonances that decay into the 
T(15) 7r + 7r~ final state, the distribution of M^p vs. 
the mass difference AM = M MAl7rw — M w is examined 
(see Fig. [T|) for the selected data sample. The cluster 
of events in the parallelogram centered at (AM, M^) 
= (1.12, 9.46) GeV/c 2 is from the transition T(45) -► 
T(lS)n + n~ . The other clusters are due to the decays 
T(2S*) -> T(lS*)7r+7r- and T(3S*) -> T(lS*)7r+7r-, where 
the T(2S) and T(3S) are produced predominantly by 
radiative return i.e. e+e" — > Y(mS)"f. 

The rightmost cluster in Fig. 1 contains events from 
the process T(45) — > T(1S)tt + it~ . The dominant back- 
ground processes, e + e~ — > /i + /i~7 (7 — > e + e~), e + e~ — > 
H + /i~H + H~, and e + e~ — > /! + /i~7r + 7r~ , accumulate at 
the kinematic boundary indicated by the diagonal line in 
Fig. 1. To capture the T(45) -> T(lS)ir + ir- signal as 
well as estimate this background more reliably, we fit the 
distribution of AM for events within the parallelogram of 
Fig. 1, whose boundaries correspond to iM^^n — y/s\ < 
60MeV/c 2 and950MeV/c 2 < AM < 1380McV/c 2 . The 
fit, shown in Fig. 2, includes a Gaussian for the signal 
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events. The ±60-MeV high horizontal shaded band is cen- 
tered on the nominal T(1S) mass. The clusters on the 
lower left correspond to T(25) — > T(15')7r + 7r _ and T(3S') — > 
T(lS)ir + ir~ transitions. The diagonal line indicates the kine- 
matic boundary A-/^ M7r7r = y/s. The parallelogram straddling 
this line defines the fitting region for T(4S) — > T(15)7r + 7r _ 
candidates. 
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FIG. 2: The fit to the AM distribution for events within 
the parallelogram of Fig. 1 using a Gaussian for the signal 
and a second-order polynomial for the background (dotted 
line). The solid curve shows the sum of the Gaussian and the 
polynomial function. 



and a quadratic function for the background. The fitted 
Gaussian is centered at (1118.7 ± 1.2) MeV/c 2 , which is 
in good agreement with the nominal 7tjy(4S) ~ m x(is) 
mass difference, and has a width of (8.1 ± 1.0)MeV/c 2 , 
which is consistent with the detector's estimated AM res- 
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FIG. 3: Background-subtracted and efficiency-corrected dis- 
tribution of 7r + 7r~ invariant mass (M,^) for events within the 
signal subregion of the parallelogram in Fig. 1. The solid 
curve shows the Mm? distribution predicted by the models of 
Ref. 0. 



olution. The signal yield in the interval 1105MeV/c 2 < 
AM < 1135 MeV/c 2 , determined as the difference of the 
number of events and the fitted background in this inter- 
val, is Nov — 113.7 ± 16.3, with a statistical significance 
of 11.2cr, corresponding to — 21n(£o/£max) = 135.0 with 
three fit parameters (mass, width, and yield). Here, Co 
and £ m ax are the likelihood values returned by the fit 
with the signal yield fixed at zero and its best fit value, 
respectively. 

Additional information can be obtained from the 
study of the ir + Tt~ system. Background-subtracted and 
efficiency-corrected distribution of 7r + 7r _ invariant mass 
(M^tt) is shown in Fig. 3 for events within the sig- 



nal subregion (1105MeV/c 2 
of the parallelogram in Fig. 



AM < 1135MeV/c 2 ) 
(The background is 



estimated from the sideband subregion 950MeV/c < 
AM < 1075MeV/c 2 and 1175MeV/c 2 < 



AM < 
with 
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14] 



1350McV/c 2 ). The EvtGcn event generator 
a matrix element that accounts for particle spins 
used to produce T(4S) — ► T(lS)7r + 7r~ — > fi + fi~ n + n~ 
events that are then passed through the detector simula- 
tion [l5| and reconstruction programs. 

The M,^ distribution in Fig. 3 can be described using 
the shape predicted by the models of Ref. [14| , in which 



suppression of small 7r + 7r invariant masses follows from 
partial conservation of axial current. The goodness of fit, 
for 10 degrees of freedom (NDF), is xV NDF = 0.35. 
The branching fraction for the T(45) — ► T(lS)ir + ir~ 



decay is determined from B(T(AS) — 
N ov /(Nr(4S) ■ e ■ B(T(LS) /i+jO), 



T(15)7r+7r-) = 
where iV PV is the 



TABLE I: Total number of T(4S r ) (N r{4S) ), signal yield (N ev ), 
reconstruction efficiency (e), and branching fraction (B) for 
the T(45) -> T(IS)tt+tv- decay. 



Data sample 




iVcv 


e, % 


B, 1CT 4 


I 


534.6±7.0 


52.2±10.7 


4.5 


0.86±0.18 


II 


122.1±1.4 


61.3±12.1 


25.1 


0.84±0.17 



of T(4S) events produced, e is the signal detection effi- 
ciency (calculated separately for samples I and II), and 
B(T(IS) -> fi+^r) = (2.48 ± 0.05)% is the PDG- 
tabulated branching fraction for the daughter decay. The 
efficiencies are calculated from Monte Carlo simulations. 
For the hadronic-cvent simulation in sample I, we apply 
a correction to E suia / yfs, one of the variables used to se- 
lect hadronic events, so that this distribution agrees with 
that of the data. This correction is also applied to sam- 
ple II, where it changes the efficiency by a few per cent. 
The results are given in Table I. 

The systematic error in the reconstruction efficiency 
due to the last correction is 8% for sample I and es- 
sentially zero for sample II. The systematic uncertainty 
in the reconstruction efficiency due to lack of knowl- 
edge of the T(45) -> T(lS>+7r- -> /i+^-Tr+Tr" de- 
cay matrix element is estimated by comparing the pa- 
rameterization of the M„ distribution in the models 
of Ref. " 



extracted signal yield, N- 



T(4S) 



is the estimated number 



14j and in a phase space model. We estimate 
this systematic uncertainty as half of the variation in 
the efficiency, and it is equal to 2.0%(3.1%); here and 
below, the first (second) value gives the systematic un- 
certainty for sample I (sample II) . The signal yield is ex- 
tracted by an unbinned extended maximum likelihood fit 
to the AM distribution for events in the parallelogram 
(Fig.l) using a Gaussian for the signal and a second- 
order polynomial for the background. The signal yield in 
the signal interval 1105 McV/c 2 < AM < 1135McV/c 2 
is determined as the difference of the number of events 
and the fitted background in this interval. The signal 
yield from the fitted Gaussian area has a larger statis- 
tical error. The systematic uncertainties from the dis- 
crepancies between these two evaluations of the signal 
yield are 2.2% and 1.2% for samples I and II, respec- 
tively. Other systematic uncertainties come from the 
choice of the fit range (0.3%, 2.5%), the choice of the sig- 
nal range (2%, 0.6%), the choice of the signal box width 
(2.4%, 1.9%), the change of the order of the polynomial 
function from two to one (0.5%, 1.2%), the tracking effi- 
ciency (4%, 4%), the muon identification efficiency (1.1%, 
1.1%), the pion identification efficiency (0.2%, 0.2%), 
the statistical uncertainty in the efficiency (1.0%, 0.6%), 
the uncertainty in the T(1S) — > decay branching 

fraction (2.0%, 2.0%), and the total number of T(45) 
events (1.3%, 1.1%). The total systematic uncertainty 
for each data sample is obtained by adding these contri- 
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butions in quadrature; the results are 10.3% and 6.6% 
for samples I and II. respectively. The systematic un- 
certainties from the tracking efficiency, T(45) counting, 
B(T(1S) — > fi + fJ.~), and muon and pion identification ef- 
ficiencies are treated as fully correlated systematic errors 
for samples I and II. Other uncertainties are treated as 
uncorrelated errors. First the weighted average of the 
uncorrelated uncertainties ((J^yg 0r ) is evaluated. The to- 
tal systematic uncertainty is obtained by adding {(J^y™ r ) 
and remaining correlated uncertainties in quadrature. 
The measured weighted product branching fraction is 

B(T(4S) -> T(lS>+7r-) x B(T(1S) -> n+n~) 
= (2.11 ± 0.30(stat.) ± 0.14(sys.)) x 10~ 6 . 
The branching fraction is 

B(T(4S) -> T(1S)tt+tt-) 

= (0.85±0.12(stat.)±0.06(sys.)) xlO" 4 . 

We also extract the partial decay width for the 
T(45) — > T(lS')7r + 7r~ transition using the world-average 
value of the total width [l[ , and obtain 

r(T(45) -► T(1S)tt+tt-) 

= (1.75 ± 0.25(stat.) ± 0.24(sys.)) keV. 

The measured values of B(?(4S) -> T(1S>+7T-) and 
r(T(45) — > T(15)7r + 7r~) supersede our previous re- 
sults with improved accuracy. The new Belle results 
are compatible with those of BABAR Q. 

To summarize, a study of transitions between T states 
with the emission of charged pions has been performed at 
Belle. The peak at AM = (1118.7 ± 1.2)MeV/c 2 is in- 
terpreted as a signal for the decay T(45) — ► T(1S)it + it~ 
with a subsequent T(1S) — ► /i + pT transition. The 
branching fraction S(T (45) — > T(lS')7r + 7r") and the par- 
tial decay width T(T(4S) —> T(lS)ir + ir~) are measured. 
We have not studied the T(45) — > T(25)7r + 7r~ decay be- 
cause criteria applied to the raw experimental data make 
our sensitivity to this decay limited. 
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